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Abstract: We study the phase diagram of SU(2) Yang-Mills theory with one adjoint Weyl 

i - 1 fermion on M 3 x S 1 as a function of the fermion mass m and the compactification scale L. 

This theory reduces to thermal pure gauge theory as m — > oo and to circle-compactified (non- 
thermal) supersymmetric gluodynamics in the limit m —> 0. In the m-L plane, there is a line 
of center-symmetry changing phase transitions. In the limit m — > oo, this transition takes 

i— i place at L c = 1/T C , where T c is the critical temperature of the deconfinement transition in 

pure Yang-Mills theory. We show that near m = 0, the critical compactification scale L c can 
J> be computed using semi-classical methods and that the transition is of second order. This 

suggests that the deconfining phase transition in pure Yang-Mills theory is continuously con- 

{S) nected to a transition that can be studied at weak coupling. The center-symmetry changing 

phase transition arises from the competition of perturbative contributions and monopole- 
instantons that destabilize the center, and topological molecules (neutral bions) that stabilize 
the center. The contribution of molecules can be computed using supersymmetry in the limit 

H m = 0, and via the Bogomolnyi-Zinn-Justin (BZJ) prescription in non-supersymmetric gauge 

^ theory. Finally, we also give a detailed discussion of an issue that has not received proper 

k^l attention in the context of N=l theories — the non-cancellation of nonzero-mode determi- 

$_i nants around supersymmetric BPS and KK monopole-instanton backgrounds on R x§ . We 

explain why the non-cancellation is required for consistency with holomorphy and supersym- 
metry and perform an explicit calculation of the one-loop determinant ratio. 
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1. Introduction 

Consider a quantum mechanical system with a potential with multiple degenerate minima. 
The ground state energy (as well as the energies of higher eigenstates) has a weak coupling 
expansion of the form 

E(g) = E pcTt . + Konpcrt. =E [1 + 0(g)] + e- 1 ' 9 [1 + 0(g)] + 0{e~ 2 / 9 ) . (1.1) 

Since e -1 / 9 has an essential singularity at g = it is impossible to express this contribution as 
a perturbative series in g, and hence this term 1 is intrinsically non-perturbative. Some of the 

lr The exponentially small terms may also be multiplied by additional negative powers and logarithms of g. 
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most interesting phenomena in quantum mechanics — tunneling, the absence of spontaneous 
symmetry breaking, the formation of energy-bands in periodic potentials — are due to e~ l / g 
effects. Although the leading term and the exponentially small contributions in (1.1) are 
intertwined in a deep way, as typical inaccuracies of perturbation theory are expressed in 
terms of functions with essential singularities as above, there is a sense in which (1.1) should 
be seen as a double expansion, a perturbative expansion in g and a non-perturbative expansion 
in e' 1 ' 9 . 

In this paper, we will use this double expansion to study the phase diagram of an asymp- 
totically free gauge theory with strong coupling scale A on I 3 x S 1 . In a theory without 
fermions the compactification scale on the S 1 circle can always be given a thermal inter- 
pretation. At small S , of size L <C A~ x , it is well-known that such theories are amenable 
to a perturbative treatment. A less widely appreciated fact is that, if certain conditions 
are satisfied, such theories are also amenable to non-perturbative semi-classical studies. Let 
$7 = P exp [i Jg! ^4^x4] denote the gauge holonomy (or Wilson line) in the compact direction, 
which, classically, is a "flat direction" . We expect that quantum effects will induce a potential 
for the holonomy Q of the form: 

V(fl) = V pcrt .(n) + Konpcrt.(^) , (1.2) 

where Vp Cr t. is the contribution of the perturbative loop-expansion in g 2 and V nonpev t.(ty is 
a non-perturbative expansion, presumably containing terms of the form e~ c ^ 9 . The pertur- 
bative term V per t. was initially computed in [1], and the calculation was extended to higher 
order in [2-4]. Although the perturbative potential V^ert.(^) is by now part of the standard 
books of thermal field theory, T-nonpert. (^) has not received as much attention. 

The perturbative calculation of the effective potential for the Wilson line in pure SU (N) 
Yang-Mills theory on M 3 x S 1 with small L = f3 gives [1]: 

2 00 1 

W«) = -Z2Q4 E nKf(l + <>tf)), (1-3) 

" 71=1 

leading to the conclusion that at small /3 the theory is in a deconfined phase, with broken 
center-symmetry (jjtrft) = 1. If one thinks in terms of eigenvalues of Q, the potential (1.3) 
generates an attraction among the eigenvalues. In other words, the effective mass-squared for 
the Wilson line is negative. 

Based on numerical simulations on the lattice we know that the deconfinement transition 
in pure Yang-Mills theory takes place at a temperature of order A: Td = aA where a is a pure 
number of order one. At one-loop order in perturbation theory, (1.3) shows that the center- 
symmetry is broken. Higher order corrections do not alter this conclusion; there is no effect at 
any order in perturbation theory that competes with center symmetry breaking. Hence, the 
phase transition must be induced by Vnonpert.(^)- Disregarding such non-perturbative effects, 
one would conclude that one cannot explore the transition as the temperature is lowered, 
from the deconfined to the confined phase, using weak coupling techniques. 

In this work, we propose a strategy to analytically study the center-symmetry changing 
phase transition in four dimensional gauge theories based on an observation discussed in [5]. 
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Figure 1: The thermal deconfmemcnt phase transition in pure Yang Mills (YM) theory can be 
accessed through a non-thermal (quantum) phase transition in supersymmetric Yang Mills (SYM) 
theory deformed by a gluino mass term. In the massless limit, the supersymmetric theory does not 
have a phase transition. The phase transition at small-™ is analytically calculable and, by decoupling, 
it is connected to thermal deconfinement phase transition in pure YM theory. 

The main idea, schematically shown in Figure 1, is as follows: It is well-known that N = 1 
SYM with periodic boundary conditions for fermions does not have a phase transition as 
a function of radius. In fact, for a supersymmetric gauge theory with Hamiltonian H and 
fermion number operator F, 

Z SYM (L) = tr [e~ LH (-l) F ] (1.4) 

is the supersymmetric (Witten) index and is independent of radius. In softly broken su- 
persymmetric theory, however, this quantity does not have an interpretation as an index. 
Consider adding a small mass for the fermion in N = 1 SYM. Eqn. (1.4) is still well-defined, 
and can be interpreted as a twisted partition function. The twisted partition function is a 
signed sum over the states in the bosonic and fermionic Hilbert spaces, Hb and Hp, according 
to the Z2 = ( — l) F grading, 

Z SYM (L,m) = Z B -Z T = e~ LEn ~ E e ~ LEn ■ ( L5 ) 

This is different from the ordinary partition function, Z SYM (/3, m) = Zg + Zjr by the over-all 
sign of the contribution of fermionic states. 

The twisted partition function, despite being a non-thermal quantity for general values 
of the fermion mass m, is immensely useful as a tool that continuously connects the thermal 
phase transition in pure Yang Mills theory with a semi-classically calculable transition on 
x Sp. A similar continuity argument at finite baryon density was made in [6]. For 
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(1.5) should be viewed as probing the phase structure of the theory as a function of radius 
L (which does not generally have an interpretation as inverse temperature). As emphasized, 
the twisted partition function is manifestly non-thermal. Yet, it can be used to study aspects 
of a genuine (thermal) deconfinement phase transition in certain limits. This is due to the 
the following decoupling argument. If the mass of the fermion is infinite, or much larger than 
the strong scale of N = 1 SYM, Z(L, m) reduces to the ordinary thermal partition function 
of pure Yang-Mills theory: 



?SYM 



(L,m) 



Z YM (/3) = tr[e- pH ] , P = L. (1.6) 



In this limit, because the heavy fermion decouples, we may identify the circumference L with 
the inverse temperature /3. For a heavy fermion, the choice of the boundary condition is 
immaterial. 

In this work, we will show that the center-symmetry changing phase transition at small m 
can be computed semi-classically. 2 In this limit the transition takes place at small L, as shown 
in Figure 1. The physics of the transition is quite interesting. It is based on the competition 
between topological molecules, called "neutral bions" or "center-stabilizing bions", and semi- 
classical monopole-instanton effects, as well as perturbative effects. We will argue that these 
effects are also present at large m, in the pure gauge theory, but that in this limit the effect 
cannot be reliably computed using semi-classical methods. 

2. Mass deformation of J\f = 1 super- Yang-Mills on S 1 x IR 3 
2.1 Perturbation theory 

Classical vacua of the theory on M 3 x S 1 are labeled by the expectation value of the Wilson 
line 



Q, = exp 



A^dxi 



(2.1) 



When LA <C 1, non-zero frequency Kaluza-Klein modes are weakly coupled and may be 
integrated out perturbatively. If we consider periodic boundary conditions for both the gauge 
fields and the adjoint Weyl fermions, A^L) = A^(0) and X(L) = +A(0), the one-loop effective 
potential for the Wilson line is [5,7]: 



7T 2 L 4 

71=1 



1+2 ( n Lm) 2 K2{nLm) 



n 4 



(2.2) 



2 In order to be precise, we note that the small-m, small-!/ calculability of the transition applies outside 
of a finite strip around the phase transition line in Figure 1. As usual in second order phase transitions, 
fluctuations become strong near the critical point and the renormalization group equations describing the 
critical theory — in the case at hand, in the 3d Ising universality class — are nonperturbative. However, in the 
weak coupling small-m, small- L regime, the width of the strongly fluctuating critical region is small, controlled 
by powers of the small parameter LA, and the critical values of the mass m c at fixed L (or critical size L c at 
fixed m) can be reliably determined up to small corrections. 
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Here m is the fermion mass and K 2(2:) is the modified Bessel function of the second kind, 
with asymptotic behavior 



As the mass m — > 00, the fermions decouple regardless of their boundary conditions, and the 
effective potential (2.2) reduces to the pure gauge result given in (1.3), with the identification 

L = f3: 

V^in,m]\ m ^ = V™[n](l + 0(e- Lm )) ■ (2.4) 

In the opposite limit of massless fermion, the one-loop potential vanishes: 

y p s e y t M [O,m = 0] = 0. (2.5) 

In fact, because of supersymmetry and the fact that (perturbatively) the theory possesses a 
moduli space of vacua, the m = SYM theory does not generate a potential for the Wilson 
line to any order in perturbation theory. At a typical point on the moduli space, a nonzero 
background Wilson line is turned on: 

e iA6>/2 

e" 



n = ( _-iA* /2 ) > (2-6) 



and the SU(2) gauge group abelianizes down to U (1). Here, AO is the separation between the 
eigenvalues of the Wilson line. Non-perturbatively, the moduli space is lifted due to effects 
leading to AO = tt, i.e., to a center-symmetric holonomy. This effect will be described in two 
complementary ways in the next Section. 

We now turn on a small mass corresponding to soft supersymmetry breaking, m <C A. In 
the small LA < 1 regime, we also have mL <C 1. In this case, using the small- z asymptote of 
(2.3), we observe that the leading term at O(m ) cancels and an 0(m 2 ) potential is induced. 
The effective potential becomes (up to 0(m ) corrections): 

(2.7) 

71=1 X ' 

where B^{x) = x 2 — x + g is the second Bernoulli polynomial (the last equality above is valid 
when AO G [0,2vr]). 

Within the domain of validity of the perturbative analysis, (2.7) shows that Wilson lines 
with all winding numbers are unstable when the fermion mass is non-zero, despite the use of 
periodic boundary condition for fermions. Consequently, the rjf = 1 theory at any non-zero 
mass m and sufficiently small L will have completely broken center symmetry. On the other 
hand, the fact that the one loop potential is small, of 0(m 2 ), implies that exponentially small 
semi-classical effects can compete with the perturbative potential. 
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2.2 Non-perturbative effects at m = via supersymmetry 

In the classical background (2.6), and at weak coupling, the Wilson line (2.1) behaves as 
an adjoint Higgs field. The theory at short distances is described by non-abelian SU(2), 
and at long distances, it is described by an abelian U(l) subgroup. Using abelian duality, 
e fj,u\d\(J = ^rF^y for /x, v = 1, 2, 3, we may map the gauge field to a spin-zero dual-photon 
a. It is also useful to define the exponent of gauge holonomy (2.6): 

b=%M . (2.8) 

t 

The kinetic terms of the fields a and b are 

In terms of superfields this corresponds to a Kahler potential: 

K = A B ' B • < 2 - 10 > 

for the chiral superfield B, whose lowest component 3 is b — ia; the fermionic component is the 
component of the gluino field A which remains massless along the Coulomb branch (2.6). The 
effective Lagrangian following from (2.10) gives the long-distance perturbative description of 
the theory on K 3 x S 1 - essentially a massless free-field theory. 

The non-perturbative dynamics of the theory is quite rich: due to the compact topology 
of the "adjoint Higgs" (2.1), there are two types of elementary monopole-instantons, M\ 
and M 2 - These are sometimes called 3d instanton and twisted-instanton, or BPS-monopole- 
instanton and KK-monopole-instanton, see [8,9]. The Nye-Singer index theorem implies two 
fermionic zero- modes for each [10, 11]. The 4d BPST instanton (in the long distance regime) 
can be viewed as a composite of these two. The monopole-instantons give rise to 't Hooft 
vertices, or amplitudes, of the form: 

Mi = e- b+ia XX, M 2 = r]e +b ~ ia XX, 

Mi = e - b - ia XX, M 2 = ve +b+ia XX, (2.11) 

a 2 

where rj = e = e~ 2So is the 4d instanton amplitude (we set the topological theta angle 
to zero). Since (2.11) carry just two zero modes, they generate a superpotential, given by: 

M% V L , R p>\ 2Ml v L „ , / 4vr 2 \ , 
W R3xS i = — ^ (e" B + 7?e B ) = e~ So cosh I B - J , (2.12) 



3 We note that the relation between B and b is, in fact, nonlinear, see eqn. (A. 35), due to the perturbative 
corrections to the moduli space metric along the Coulomb branch. These are also reflected in the non- 
cancellation of the one-loop fermion and boson nonzero mode determinants around the BPS and KK monopole- 
instantons, see Appendix A for a detailed discussion. As these subtleties represent subleading corrections to 
the Kahler metric (2.10), we ignore them in the main text. 
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where the coupling is normalized at the cutoff scale Mpy] details of the instanton calculation 
leading to (2.12) can be found in [12]. Thus, the infrared Lagrangian is given by (2.10) and 
(2.12): 

£= [ d*6 K + ( f d 2 9W + h.cA , (2.13) 



and the scalar potential can be easily found: 

t,„ \ r _i dW 2 64 7 r 2 M^ / L 3 e- 2So / , / /, 4vr 2 ' 
V(b,a) = K BtB 



dB 



'^-^ (cosh ( 2 (b-^-Jj -cos 2.), (2.14) 



where i^BtB i s the mixed second derivative of the Kahler potential. Furthermore, it is con- 
venient to introduce, instead of b of eqn. (2.8), the shifted field: 

Finally, to rewrite the potential in terms of the strong-coupling scale A and the S 1 size L, we 
use the relations: 

^«-*-A>. ^"Slog* (2.16) 
3 g 2 AL 

These relations express the fact that the scale for the coupling of the effective theory is set by 

the compactification scale L. We use the one- loop running coupling constant in the prefactors 

and two-loop running in the exponent of the instanton amplitude. 

Thus, using (2.15) and (2.16), we obtain the final expression for the scalar potential 

(2.14): 

V hion (b', cr) = 48L 3 A 6 log — (cosh 2b' - cos 2a) . (2.17) 

Note that if the superpotential W is determined by BPS and KK monopole-instantons then 
the scalar potential is governed by (correlated) monopole-anti-monopole pairs. We have 
therefore denoted the scalar potential by Vbi on - We will make this relationship more explicit 
in the following Section. 

The potential Vbion for the ft'-field, ~ cosh(26'), generates a non-perturbative repulsive 
interaction between the eigenvalues of the Wilson line around S 1 : it is minimized at (b') = 0, 
which, from (2.15), corresponds to maximally separated eigenvalues, {AO) = tt, or in terms 
of the Wilson line: 

(n) = (e iA ± L ) = _,A (trfi)=0, (2.18) 



up to gauge rotations. This is the center-symmetric vacuum of the theory on R 3 x S . On 
the other hand, the Vbion potential for the a field has two minima, located at 

(<7)={0,7T} 5 (2.19) 

associated with discrete chiral symmetry breaking. Evidently, in the effective lagrangian, the 
mass gap for gauge fluctuations a is generated by the operator e ±2j<J , and for the spin-zero 
scalar (the fluctuation of gauge holonomy b') it is generated by e ±2b '. 
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Expanding the action around the center-symmetric gauge holonomy and using component 
notation, we find the effective lagrangian: 

1 9(L) 2 Ua u „ 2 , , 21 , , L 



C= ^j^Z Wf + idiaf] +i^Xa i d i X + ae-^W y(e~ b ' +io ' + e +b '~ icr j XX + c.c. 

— „ 

+ p e — F ^\e- 2b ' + e 2h '-e- 2l °-e 2i °] . (2.20) 

We expressed the Lagrangian in component notation to elucidate the physical origin of the 
various terms 4 . 

2.3 Non-perturbative effects at m = via topological molecules 

We now provide a derivation of the bosonic potential without the use of supersymmetry. This 
formalism will apply to both supersymmetric and non-supersymmetric theories. As discussed 
in Refs. [13,14], in general gauge theories with massless adjoint fermions (i.e., not only in 
the supersymmetric single massless Weyl flavor case rif = 1), the 3d instanton and twisted 
instanton do not generate a mass gap for the gauge fluctuations because of their fermionic 
zero mode structure. The zero modes are determined by the index theorem [10, 11], and 
the corresponding instanton amplitudes have the form given in (2.11) or the generalization 
thereof for the nf > 1 case. 

Unlike the superpotential, which arises due to monopole-instantons with exactly two 
zero modes, the associated bosonic potential must be induced by topological molecules which 
do not have any fermionic zero modes. The bosonic potential is generated by correlated 
monopole-anti-monopole pairs. To second order in the semi-classical expansion, the possibil- 
ities are the following: [.M1.M1], [M2M2], [M1M2], [M2M.1]. These objects can be viewed 
as composites of (2.11). The magnetic and topological charges and the amplitudes associated 
with these instanton-anti-instanton events are: 



composite 


(Qm, Qt) 


amplitude 


[M x Mx] 


(0,0) 


e~ 2b ' 


[M2M2] 


(0,0) 


e +2b ' 


[M{M 2 ] 


(+2,0) 


e +2la 


[M2M1] 


(-2,0) 


e" 2lCT 



(2.21) 

The action and interaction due to massless boson exchange of two fundamental monopole 
instantons with (bosonic) amplitudes e n & fc '+ m ™°" an d e n 'f> b ' +m ' m(7 is (here rib,n m = ±1): 

S(r) = 2S + S int = + (-n b n' b + n m n' m ) (2.22) 

gZ gZj, 

4 In order not to clutter notation, we kept only the exponential dependence of the coupling g 2 (L) in the 
non-perturbative terms in (2.20), i.e., we absorbed the numerical coefficients and the log -r dependence in 
the prefactors a and /3; these can be recovered from (2.f2), (2.17). 
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where r = |r| is the separation between two instanton events. The interaction proportional 
to n^n! h is due to the exchange of the massless A4 modulus and the one proportional to 
n m n' m is due to the exchange of the dual photon ("magneto static"). The interaction term 
is repulsive for (— n^nl + n m n' m ) > 0, attractive for {—n\,n' b + n m n' m ) < and zero otherwise. 
The fermion-induced interactions will be considered below. 



2.3.1 Magnetic bions, quasi-zero modes, and the mass gap for the dual photon 

All the topological molecules contributing to the bosonic potential have vanishing topological 
charge, i.e., they are indistinguishable from the perturbative vacuum in that sense. However, 
the [.M1.M2] (and its anti-molecule) events carry two units of magnetic charge. The prefactor 
of the amplitude can be found as follows. 

The Mi instanton has four bosonic zero modes. Three of these are the positions xeR 3 
and one is related to the internal U(l) symmetry. Note that there is no size modulus associated 
with monopole-instantons, unlike the 4d BPST instanton. This is one of the reasons that we 
can do reliable semi- classical analysis. Let an Mi be located at y. Because of the interaction 
which depends on the separation of the two events, r = x — y, the relative coordinate is no 
longer an exact zero mode, while the "center of mass" position R = (x + y)/2 of the two 
events is an exact zero mode. The relative coordinate is now a quasi-zero mode. This is to 
say that the operator of the quadratic fluctuations in the background of Mi and M2 has, in 
its spectrum, an exact zero mode, a low lying quasi-zero mode and parametrically separated 
Gaussian fluctuations. The latter modes can be trivially integrated out. The zero and quasi- 
zero modes are particularly important. In particular, the integrals over the quasi-zero modes 
need to be done exactly. 

The magnetic bion amplitude associated with a bion located at R can be found by 
integrating over the quasi-zero mode exactly. Below, we write the expression for the rif flavor 
theory for later convenience. The amplitude associated with an [.M1.M2] composite is: 

[M1M2] ~ Ae- 2So e 2ia , (2.23) 

where 

A = [ d 3 re-( 2x ^ +4 "' losr ) = 47r/(A,ra/) , A ee /- . (2.24) 



The meaning of the terms in the exponent is as follows: 2 x accounts for the repulsion 
due to exchange of a and 6-scalars, and 4njlogr is the attraction due to fermion zero mode 
exchange. Consequently, there is a single saddle-point in the quasi-zero mode integral, given 
(for rtf = 1) by: 

r b = ^ , (2.25) 

g 2 

which can be interpreted as the magnetic bion size. The bion size is much larger than 
monopole-size, but much smaller than (uncorrelated) inter-monopole separation. Conse- 
quently, a representation of the partition function as a dilute gas of magnetic bions is justified. 
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The integral in (2.24) is given by: 

J(A,n/) = J dre-^+^n s -2)\o g r) = X in f~ 3 T{4n f - 3), (2.26) 

where nf = 1 for SYM. The way to check that this is self-consistent is as follows. The 
interaction term for the magnetic bion in (2.22) must be parametrically smaller than the 
leading action in order for the relative coordinate to deserve the name quasi-zero mode. 
Just saying that 25o S> S- mt , or equivalently, r 3> ^, is not sufficiently good, because this 
does not preclude the O(l) changes in the combined action, 2Sq. It must be such that 
|Sint| ~ 0(g 2 ) x (2,So) <C (2So). Indeed, the action of the magnetic bion configuration 
associated with r = r^, given in (2.22), takes the form 

S(r h ) = 2S (1 + 0{g 2 )) (2.27) 

Since the interaction changes the action only by parametrically small 0(g 2 ) effects, the mag- 
netic bion topological molecule can be viewed as a quasi-solution. 

The magnetic bion molecules described in this Section are responsible for the generation 
of a mass gap for the dual photon (the potential for a in (2.20)) in SYM on M 3 x S 1 , and they 
generate the confining string tension. In the following Section we will study a second type 
of topological molecule, which is more subtle to identify, but plays an important role in the 
center-symmetry realization. 

2.3.2 Neutral (center-stabilizing) bions and the BZJ prescription 

Consider now the other possible composite from the list (2.21), the [.M composite which 
carries no magnetic and topological charge (the [.M2.M2] is treated similarly). Here, the 
integral over the quasi-zero mode is, naively: 

[M1M1] ~ Aaiv C e- 25o e ±2b ', (2.28) 

where: 

AaivcG? 2 ) = / d 3 r e -(- 2x ^ +4 ^ logr ) =47r/(A,n/) , 

I(X,n f ) = J dre-(-^ +(4n '- 2)logr ) . (2.29) 

Now, the interactions between constituents due to a and b exchange are both attractive, while 
the fermion zero mode induced attraction is not altered (it remains attractive). The integral, 
(2.29), as it stands, is dominated by the small r regime, where not only (2.29) is incorrect, it 
is also hard to make sense of constituents as the interaction becomes large. This is in sharp 
contrast with the magnetic bion (2.26) [14,15]. 

At first sight, this may seem to prevent us from computing the contribution from these 
pairs, but this is not actually the case. We will reach a satisfactory resolution of the problem, 
via the Bogomolnyi-Zinn- Justin prescription: The integrals over the quasi-zero modes of at- 
tractive instanton-anti-instanton molecules can be calculated in a manner initially described 
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by Bogomolnyi [16] in the context of quantum mechanics. The relation between this pre- 
scription and the Borel procedure was pointed out in the same context by Zinn- Justin [17], 
see also [18]. The prescription is to analytically continue the coupling g 2 in the instanton- 
anti-instanton interaction to — g 2 . This turns the attractive Coulomb force into a repulsive 
one. Then we calculate the resulting integral exactly, without any Gaussian approximations. 
Finally, we analytically continue the final result back to positive g 2 . 

Recall that, very often for a non-Borel summable series, when g 2 is continued to — g 2 , 
the series become Borel summable. However, one needs to continue back to positive g 2 . 
Depending on the path that one takes the coupling to the positive g 2 , the Borel sum typ- 
ically produces an ambiguous (non-perturbative) imaginary part. This is a manifestation 
of non-Borel summability. In the quantum mechanical examples that Refs. [16] and [17] 
studied, the ambiguity in the Borel sum is canceled by the ambiguity associated with the 
attractive instanton-anti-instanton molecule. The prescription for the topological molecules 
may be viewed as consistently extending the Borel prescription for perturbative sums to 
non-perturbative sectors with vanishing quantum numbers. For a fuller discussion of these 
phenomena in field theory, see [19]. 

This prescription will give an overall phase between the magnetic bion amplitude and 
center-stabilizing bion amplitude. This phase difference is physical and crucial for our con- 
siderations. Following the prescription, when we modify /(A, n/) — > I(—X,nf) (recall that 
A = g 2 /8itL) the Coulomb- interaction becomes repulsive and we can evaluate the integral 
over the quasi-zero mode. In fact, it is equal to the integral I(\,rif) for the magnetic bion 
computed in the previous Section. Next, we substitute g 2 — > —g 2 (A — > —A) to obtain the 
final result for the center-stabilizing bion amplitude. To summarize, the generalization of the 
BZJ prescription to field theory results in the chain: 

I(X,n f ) -+ I(-X,n f ) = I(X,n f ) -+ /(-A, n f ) = (-A) 3 " 4 "' T(4n f - 3) = -/(A, n f ) . 

(2.30) 

The last equality is only valid for integer n / and gives an overall sign of the center-stabilizing 
bion amplitude opposite that for the magnetic bion. We note that this line of reasoning has 
a close parallel in supersymmetric quantum mechanics [21]. 

The importance of the relative sign between the magnetic bion amplitude and center- 
stabilizing bion amplitude is worth noting, as it is a physical consequence of our prescription. 
As a result, we obtain for their combined contribution: 

V(b,a) ~ r/ cosh 26' - r/ cos 2a = e~ 2S ° [(1 + 2b' 2 + ...)- (1 - 2a 2 + .. .)] 

= 2e- 2So (b' 2 + a 2 ) , (2.31) 

the same result that we obtained earlier by using holomorphy. The crucial point here is 
the cancellation of the "cosmological constant" term in the potential. Recall that in a su- 
persymmetric theory with unbroken supersymmetry, the expectation value of Hamiltonian is 
positive semi-definite and (\l/ n |.ff|vl/ n ) > 0, and that the bound is saturated for the ground 
state (^ol^l^o) = 0. If the relative sign was not present, the ground state energy would not 
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vanish, implying a breakdown of supersymmetry. Equally importantly, the absence of the 
relative sign between the [-M1-M1] and [.M1.M2] amplitudes would lead to the presence of a 
relative sign between the mass term for the two scalars of the form (b' 2 — a 2 ), signaling an 
instability. Clearly, neither is the case. 

2.4 Center-stabilizing vs. center-breaking effects in softly broken SYM 

As explained in Section 2.1 there is no perturbative contribution to the Wilson line effective 
potential for m = 0, but there is a non-perturbatively induced potential which ensures unbro- 
ken center symmetry in the supersymmetric theory on I 3 x S 1 . This potential, as explained 
above, is due to center-stabilizing bions. 

We now turn on a small but non-zero m. At small L, one expects a competition between 
the one-loop 0{m 2 ) potential for the Wilson line and the non-perturbatively induced super- 
potential, leading to non-uniformity in the m — > and L — > limits. Taking m — > first, 
the theory lands on the center-symmetric phase. If, instead, the L — > limit is taken first, 
the theory lands on the center-broken phase. The transition line separating center-symmetric 
and center-broken phases must emerge from the L = m = corner of the phase diagram, as 
illustrated for an SU(2) theory on Fig. 1. Let us now describe the center-symmetry breaking 
dynamics in some more detail. 

Adding a soft mass term for the fermions reduces the J\f = 1 supersymmetry to J\f = 
and has the effect of lifting the fermion zero modes from the instanton amplitudes (2.11). 
The mass perturbation is: 

TTl 

AC rn = -5 trAA + h.c. . (2.32) 

g 2 

The insertion of mass terms lifts the zero modes of the monopole-instanton amplitudes (2.11) 
which now contribute to the potential for a and b. The corresponding calculation is presented 
in the appendix, and the result for the monopole-instanton contribution to the scalar potential, 
to leading order in m (mL <C 1), takes the form: 

V mon _ = 24mLA 3 cos a ^log cosh b' - ho' sinh b'^j . (2.33) 

Despite the addition of a fermion mass term, the fermion-attraction mechanism giving rise to 
magnetic and center-stabilizing bions is still operative, provided the fermion mass is smaller 
than the inverse size of the bions, = from (2.25), i.e., for mL < |^ (below, we 
shall see that this condition is obeyed in the regime where we can study the competition 
between center-breaking and center-stabilizing effects). Thus, adding the bion and monopole 
non-perturbative contributions (eqns. (2.17) and (2.33), respectively) to the the perturbative 
contribution (2.7), recalling (2.15), we obtain the full scalar potential: 

Vtatel = 48L 3 A 6 log — (cosh 2b' - cos 2a) 

+ 24?nLA 3 cos a (log -L cosh b' - h' sinh If) ^ \ (b') 2 . (2.34) 

V AL 3 / 36L log 
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It is convenient to introduce dimensionless masses, compactification scale, and potential: 



~_ m f_ A r t7— -^ 3 ^total , oc .\ 

m = — , L = AL , V = — ~ ~ — . (2.35 

A 48L 6 logL^ 1 V ' 

The final result for the scalar potential of the mass-deformed SYM theory is: 

V = cosh 2b' — cos 2a H — ^— cos a (cosh b' ~ — b' sinh b' 

2L 2 V SlogL" 1 

-I^(l) 2 i^FT (6 ' )! - (2 ' 36 » 

The physics that this potential encapsulates is our main result. Before we study the 
relative importance of the various terms in (2.36), let us summarize the region of validity 
of the scalar potential. It was derived using weak-coupling semi-classical calculations at 
small L, whose validity requires that AL <C 1. The validity of the fermion-pairing bion 
mechanism further requires mL log -rj <C 1 ; the usual soft-breaking condition m <C A is then 
automatically satisfied. In other words, both dimensionless parameters rh and L from (2.35) 
are small. 

We can now use the potential (2.36) to study the symmetry realization of the theory as 
the parameters are varied: 

1. Consider the domain: 

rh 

« 1 , (2.37) 

In this domain, the bion-induced center-stabilizing term (cosh 2b 1 ) dominates over both 
the monopole and perturbative contributions, both of which favor center-symmetry 
breaking, as we show below. In the regime of small -5-, with rh > (recall that we set 
the #-angle to zero), the vacuum with (a) = ir, (b') = represents the global minimum 
of V, while the one with (a) = is only metastable. In this regime, (trQ) = 0, with 
unbroken center symmetry. To see this explicitly, we expand (2.36) around (a) = ir, 
(b') = 0, to quadratic order in the fluctuations 5a and 5b': 



V = 2\l-^ 
8L 2 



2ml 

1 + 



31ogL L 2 432 log L 



(5b') 2 



m 



+ 21 1 + 1 (Saf , (2.38) 

where we have dropped the constant V(b' = 0,a = tt). In the regime (2.37) there is a 
mass gap for the dual photon and the theory is in the confining phase with unbroken 
center-symmetry. 

2. It is clear from (2.38) that, as we depart the regime (2.37), center-symmetry becomes 
destabilized. The leading center-breaking effect, in the region where the semi-classical 
analysis is valid, is due to the non-perturbative monopole-instanton term, which tends 
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to reduce the mass of the Wilson line. The effect of the center-destabilizing perturba- 
tive contribution, given by the last term in (2.36), is suppressed both numerically and 
parametrically, by the large factor | log 3 L\. 

Thus, keeping to small and fixed, as we further decrease L, the theory leaves the con- 
fining domain (2.37). The monopole term becomes the most dominant and destabilizes 
the (&') = center-symmetric vacuum. As eqn. (2.38) shows, the center-symmetry 
destabilization is continuous. At L = L c where 



the center-symmetric vacuum gets destabilized and the two eigenvalues of the holonomy 
smoothly approach each other. 

The potential (2.36) can be used to study the physics until the semi-classical description 
breaks down (this occurs when the scale of SU(2) — > U (1) breaking governed by the eigenvalue 



the appearance of two vacua. In the center broken phase, we expect that these two vacua are 
continuously connected to the two thermal equilibrium states of pure Yang-Mills theory as 
rh — > oo. 

One crucial point here is the following. In the confined phase, the effective description of 
the dynamics is given in terms of the Wilson line O (the 6'-field) and dual photon a. On the 
other hand, a is not a well-defined notion in the "deep" deconfined phase where the SU (2) 
gauge symmetry is fully restored and the abelianization of the dynamics is lost. In other 
words, the combined potential (2.36) is strictly valid beyond L > L c (m) and for a range 
L < L c (m) provided the eigenvalues are sufficiently apart. For most of the range L < L c (m), 
the potential is solely in terms of fi, without a. 

Finally, we can try to perform a (very rough) extrapolation of our result to pure Yang- 
Mills theory and obtain an estimate of the critical temperature of the deconfining phase 

transition. In the semi-classical domain, from (2.39), we find j£ = ~ v/f^> w hich drops 
with to, but for to 3> A the result must become independent of to. Not much is known 
numerically about the decoupling scale for a Weyl fermion in the adjoint representation. In 
the case of N c = 3 QCD with three flavors of fundamental fermions it is known that relatively 
large values of the fermion mass, to > 5A, are needed in order for the phase transition to 
approach the deconfinement transition of the pure gauge theory [22]. Assuming that the 
decoupling scale for an adjoint Weyl fermion is in the range to^ gc ~ (5 — 10)A we expect 
T c ~ (0.8 — 1.3) A, broadly consistent with lattice data. 

3. Pure Yang-Mills theory 

In the previous Section, we showed that for to <C 1 the center-symmetry restoring phase tran- 
sition can be described semi-classically. In this regime the transition is driven by the com- 
petition between center-stabilizing topological molecules and center-destabilizing monopole- 
instantons. In this Section, we will show that the same mechanism also exists in the pure 




(2.39) 



difference AO times £ becomes comparable to the strong-coupling scale A). The evolution of 
the eigenvalues leads to spontaneous breaking of the Z2 center symmetry, (itrfi) = ±1, and 
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gauge theory, even though in this case the effects cannot be computed reliably. This implies 
that it is plausible that the deconfinement transition in pure gauge theory is driven by the 
same topological phenomena that operate in the small rh limit. 

3.1 Non-perturbative effects on the classical background 

In this Section, we will consider possible non-perturbative contributions to the potential for 
the Wilson line in pure Yang Mills theory. The question is whether there are terms that favor 
the center-symmetric vacuum and compete with the perturbative contributions to V(fl). We 
consider a classical background field on M 3 x S 1 : 



iA6/2 



e -iA0/2), (3-1) 

where A9 is the separation between the eigenvalues of Wilson line. In the classical background 
(3.1), and at weak coupling, the Wilson line behaves as an adjoint Higgs field breaking the 
microscopic SU(2) symmetry down to U(l) at large distances. As before there are two types 
of elementary monopole-instantons, Mi and Mi- The amplitudes associated with these 
instanton events are essentially the ones given by (2.11), but now without the fermion zero 
modes: 



Mi = e F"^- = e ~ b+la , Mi = e ^ = e 

(3 2) 

M2 = e s 7 = rje +b la , M2 = e ^ = rje +b+lcr . 



The interaction between different monopole-instantons with magnetic charge n m and scalar 
charge rib can be computed by using the two point correlator, as in (2.22): 

/ e n b b+in m cr f x \ e n b b+ih m a _ g-V(|x-y|) _ g - g Z%ty\ (-™6«6+™mn m ) ^ ,^ ^ 

This formula is true classically. At weak coupling, the b field may acquire a radiatively induced 
mass. This modifies the potential as: 

F(|x - y|) = 2| 47rL (-rw-^-yl + n m n m ) . (3.4) 
9 |x-y| 

In perturbation theory, there are three possibilities for (recall that is the mass for the 
Wilson line, m||trf2| 2 , obtained by studying small fluctuations around the classical vacuum): 
i) In supersymmetric theories with supersymmetry preserving boundary conditions, = 
to all orders in perturbation theory, ii) If > 0, as is the case for QCD(adj) with periodic 
boundary conditions for fermions [13,14], then the potential is as in (3.4). In particular, the 
6-exchange interaction is short range and has no effect on the long distance effective theory. 
in) If m b < 0, as in thermal YM, then there exists an instability of the center-symmetric 
vacuum. 

For our purposes we will consider = 0, because classically there is no mass gap for 
the holonomy fluctuations. We can now write down the effective theory for an ensemble of 
monopole-instantons. We get: 

L = [m 2 + (d t a) 2 ] - (e~ b + r,e+ b ) «** + ..., (3.5) 
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where we have not attempted to determine the overall coefficient of the effective potential, 
and we defined b = ^AO. The ellipsis denote both perturbative and other non-perturbative 
contributions. The effective potential in (3.5) arises by summing the contributions of the 
different monopole-instantons in (3.2): 

-V n - p -(A8,a) = (M!+M2 + M2 + M 2 ) 



[AO -i£(2n-A0) 



cos a = 4e "a 2 " cosh b' cos a . (3.6) 



We would like to make a number of comments regarding this potential: 

1) We note that the potential is not quite sensible as it is not manifestly periodic in AO. 
This problem can be addressed by taking into account the existence of an infinite tower of 
monopole-instantons, see below. 

2) We observe that the potential has an extremum at a = and AO = where the 
semi-classical approximation is not reliable. For the stability of the center-symmetric vacuum 
we must have a global minimum at AO = tt, and the sigma field must be non-tachyonic there. 
However, around (AO, a) = (vr,0), and also around (AO, a) = (it, it), the expansion of the 
potential in terms of small fluctuation yields: 

[ -(5b') 2 + (5a) 2 + ... for (AO, a) = (tt, 0) 
V n - p -(b',a) = 2e ^ { ' \ ■ (3.7) 

' +(5b') 2 - (5a) 2 + ... for (AO, a) = (tt, tt) 

At leading order in the semi-classical expansion, the Hessian around each of the two center- 
symmetric saddle points is negative and we conclude that monopole-instanton effects do not 
favor the center-symmetric vacuum. 

In the following we will address both of these points. We will argue, in particular, that 
topological molecules can stabilize the center. 

1) Making the potential periodic: Because the AO field is an angular variable there is an 
infinite set of monopole-instantons on E 3 x S 1 . For magnetic charge +1, there is a tower of 
monopole-instantons with topological charges QT( n w) = + n w ,n w S 7L. We observe that 
Mi and M2 are the n w = and n w = — 1 members of this tower. Following [23] we define 
the generalized fugacity of the monopole-instanton amplitude as: 



F(A0)e i 



E 



-^\Ae+2nn w \ 



e ia , (3.8) 



e » 

Clearly, F(A0 + 2ir) = F(A0) is a periodic function. Using Poisson resummation it can be 
rewritten as: 

47T / OO \ 4-7T 



F(A0) = -J2 ^ ineAd = ~ 1 + 2 E c ^ cosn eA0 , c n 

7r n e ez[H) +n 2 n V n e =l / 



which, after using cos n e AO = i|trft ne | 2 - I, can be expressed in terms of the gauge invariant 
Wilson line: 



F(fi) = Ml+f;c ne (|tr^| 2 -2)) . 

\ n e =l / 



(3.9) 
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The non-perturbatively induced potential, incorporating the entire Kaluza-Klein tower of the 
monopole-instantons takes the form: 



(3.10) 



The extremum of (3.10) is at A8 = 0, a = it. At this point center-symmetry is broken, gauge 
symmetry is restored and the abelian long distance description is invalid, exactly as for (3.6). 
Thus, the leading-order bosonic potential induced by monopole-instantons does not appear 
to stabilize center symmetry. 

2) Role of topological molecules: We showed that in J\f = 1 SYM 5 monopole-instantons 
carry fermionic zero modes and do not induce a bosonic potential. Instead, monopole- 
instantons generate a fermion mass term. A bosonic potential appears at second order in 
the semi-classical expansion through the terms induced by composites [.Mi-Mj]. The mag- 
netic bion combinations [.M1.M2] + [A^-Mi] induce a mass gap for gauge fluctuations, and 
the center stabilizing bions + [.M2.M2] stabilizes the center-symmetric vacuum. We 

may therefore ask whether topological molecules induce similar effects in pure Yang-Mills 
theory. 

Symmetry permits, at second order in the semi-classical expansion, terms of the form: 



The terms in (3.11), which can be made periodic in b by a procedure similar to the one that 
led to (3.10), can be thought as due to molecular monopole-instantons of the form: 



In Section 2.3.2, we showed that there are SUSY and softly broken SUSY theories in which we 
can reliably demonstrate that these molecules generate a center-symmetric minimum. Below 
we will demonstrate, using the BZJ prescription, that topological molecules also stabilize the 
center in pure YM theory. Clearly, in order for the second order terms in (3.11) to be more 
important than the leading-order term in (3.10), we have to push the expansion beyond the 
semi-classical domain. We will therefore not attempt to perform a quantitative calculation. 
Our goal is to show that the same mechanism that drives the center-symmetry changing 
transition in the softly broken M = 1 theory also operates in pure Yang-Mills theory. 

3.1.1 Quasi-zero modes and bion amplitudes 

In this Section, we will study the amplitude of topological molecules in the pure Yang-Mills 
theory on the classical background. The amplitudes for the molecular monopole-instanton 

5 The same is true in Af = 2 Seiberg-Witten theory softly broken to M = 1. 




[M1M2] ~ e +2i ° , 
[M1M1] ~ e~ 2b , 

[M1M2] ~ V , 
[MiMx] ~ e~ 2b+2ia , 
[M2M2} ~ v 2 e +2b - 2iCT 



[M2M1] ~ e~ 2ia , 
[M2M2] ~ V 2 e +2b , 



[M1M2] ~ 7] , 

[MiMi] ~ e~ 2b - 2ia , 
[M2M2] ~ i 1 2 e +2b+2iu 



(3.12) 
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events can be found by integrating over the relative separation quasi-zero mode between 
its constituents. This will be similar to the study in Section 2.3.2, where we performed an 
analysis for the case of nj adjoint fermions, where nt = 1 corresponds to N = 1 SYM theory. 
The result for pure Yang Mills theory can be obtained by taking the limit nt = e — > in 
(2.26), (2.29), and (2.30). 

We begin with topological molecules of the type [.M1.M2]. The integral over the quasi- 
zero mode is the same as in (2.26), except that we have to take the e — > limit. We find 

/(A, e) = A 4 - 3 r(4e - 3) = A" 3 (- ±- + ± (- log(A) + 7 - y ) + 0(e)) • (3-13) 

The divergence in the e — > limit is due to over-counting of uncorrelated monopole-instanton 
events, which are already included in the dilute monopole-instanton gas approximation. In 
theories with fermions, this long-distance divergence is cut-off by the fermion zero mode 
exchange, both in quantum mechanical examples [21] and in quantum field theories [19, 
23]. In theories without fermions, in order not to double-count, this divergence needs to be 
subtracted, see [16, 18] for a quantum mechanical example. Consequently, the prefactor of 
the magnetic bion amplitude is: 



*«-^r?y(-*(£W!M. 



where a is the coefficient of the one monopole-instanton amplitude, which is set to one in 
(3.2) and which can be restored if desired. 

Next, we consider the neutral bions [«MiA4i]. Here, the constituents interact attractively 
both due to a and 6-exchange, and we need to apply the BZJ prescription. The result is: 

, s 2na 2 ( 8tt\ 3 / , ( 9 2 \ ll\ , N 2vra 2 / 8vr\ 3 . 

In analogy with the quantum mechanics example where the ambiguity associated with non- 
Borel summability of the perturbation theory is canceled by the molecular instanton-anti- 
instanton contribution [16,17], we also expect the ambiguity associated with the non-Borel 
summability of the pure YM theory on M 3 x S 1 to be canceled by the two-fold ambiguity of 
the amplitude for neutral bions. This implies that the imaginary part in (3.14) must cancel 
by large-orders in perturbation theory, and we discard that term. The remaining term in 
(3.15) satisfies 

C2{g) = -ci(g) . (3.16) 

For a more detailed discussion of this conjecture and the available theoretical evidence in its 
favor, see [19]. 

Next, consider the [A^iA^i] and [.A^-A^] amplitudes. Now, a exchange is repulsive and 
b exchange is attractive, leading to a cancellation. This means that the self-dual monopole- 
instanton events are not correlated, and the molecular amplitude vanishes, c^,{g) = 0. The 
same conclusion is also obtained in Ref. [24]. 
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For the LM1.M2] amplitude, slightly more care is needed. The interaction again cancels 
between attractive a exchange and repulsive b exchange. This would seem to imply 04(g) = 0, 
as in the previous case. However, [.A/f 1.M2] molecules with sizes less than the compactification 
scale correspond to ordinary BPST-instantons. Taking into account small 4d instantons cor- 
responds to 04(g) ~ 77. Large 4d instantons do not exist because large uncorrelated [.M1.M2] 
molecules are already included in the instanton-monopole contribution. This implies that 
there is a finite instanton term, but since the 4d-instanton does not depend on the b and a 
field we find that the leading semi-classical instanton contribution only enters as a constant 
term in the effective action. 

Combining the terms (3.14), (3.15), and (3.16) that appear at second order in the semi- 
classical expansion we find: 

V< 2 ) (b, a) = -cir ? cos2 ( j + ci (e -26 + 77V 6 ) . (3.17) 

This is indeed the same result as in (softly broken) N = 1 supersymmetric theory. The 
combined potential of monopole-instantons and bions is 

V {1) (b,a) + V [2) {b 1 a) = -r/ 1/2 cosh b' cos a - c x r] cos 2a + c\T) cosh 2b' , (3.18) 



leading to: 



ml 



ml 



-rj 1/2 + Ac\r) 



X ' 2 + 4 Cl 



(3.19) 
(3.20) 



Thus, we find that second-order effects in the semi-classical expansion, in particular neutral 
bions, stabilize the center symmetry, without destabilizing the a mode. As in the discussion 
around (3.9) and (3.10), in order to make periodicity of the potential manifest, we may sum 
over the Kaluza-Klein tower of the neutral bions. This leads to a center-stabilizing potential 
for the Wilson line, given by Vb-tower(^) ~ Yin |trf2 n | 2 , similar to (3.9). However, as stated 
earlier, the stabilization of center-symmetry by the neutral bion induced potential requires 
strong coupling where neutral bion term can overcome the monopole-instanton induced in- 
stability as well as the perturbatively induced potential. In this regime, we cannot perform a 
quantitative calculation. It is nevertheless intriguing that the same effect that leads to center 
stabilization at weak coupling for softly broken M = 1 theory is present in pure Yang Mills 
theory as well. For an attempt to connect semi-classical neutral bion molecules to strong 
coupling effects and infrared renormalons, see the recent work [19]. 



3.1.2 Connecting monopole-instantons to dyon particles 

The neutral bion induced potential (3.17), just like the monopole-instanton induced potential, 
is not manifestly periodic in AO. This problem, as in the case of monopole-instantons, can 
be addressed by taking into account the whole Kaluza-Klein tower of neutral bion molecules. 
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In this Section, we will show that the tower of monopole-instantons can be interpreted, via 
Poisson resummation, as the contribution from dyon particles. 6 

Consider the Poisson duality relation for the tower of monopole-instantons, see (3.9): 



e a 



\A9+2nn w \ 



An e A9 



(3.21) 



The phase e tneLA4 = e me ^ 9 is the canonical coupling of a charged particle to a background 
gauge field A4, and e lu is a 3d instanton amplitude associated with magnetic charge one. 

Consider a massive particle on M 3 ' 1 with magnetic and electric charge (n m ,n e ) and mass 
M( nm ne ), called a dyon. By Gauss' law, there is a combined electric and magnetic flux coming 
out of this particle; 



(E + iB) ■ dT, = n e + i^n m = q e + iq. 



(3.22) 



where S 2 = <9M 3 is a sphere at infinity. We can analytically continue to Euclidean time 
(consider the theory on M 4 ) and then compactify one of the directions, i.e., consider the 
theory on IR 3 x S 1 . 

A dyon particle whose world-line wraps the S 1 corresponds to a finite action instanton, 
S(n m ,n e ) = ^^(n m ,n e )- The amplitude associated with the dyon is, 



-LM 



(n m ,n e ) p ine8 e +in m cr 



(3.23) 



where e l ® e and e J0m are electric and magnetic Wilson lines. The magnetic Wilson line is 
naturally interpreted in terms of the dual photon 7 a. The duality relation (3.21) therefore 
admits an interesting interpretation: its right hand side may be re-written as a sum over 
massless/light dyons with charges (n m ,n e ) = (l,n e ) whose world-lines wrap around the S 1 : 



1 



7T 



—LM(q m ,q e ) in e 9 e +iT 



(3.24) 



6 A number of authors, see, e.g., [20] and references therein, refer to the monopole-instantons Mi and AI2 
as "dyons", because they are magnetically charged and self-dual. We believe that this choice of words is not 
quite appropriate. The amplitude of a monopole-instanton has the form e _i>+1<T ~ e _A4+I<T . However, a dyon 
particle with both electric and magnetic charge couples to (Ai,a) as e iq ' Ai+iqma ', in both Minkowski and 
Euclidean space. We will see that the Kaluza-Klein tower of monopole-instantons can be written as a sum 
over dyons that exhibit the expected coupling for electrically and magnetically charged particles. 

7 We can see this starting from abelian duality in 4d, then compactify the theory on R 3 x S 1 , and finally 
match the result to 3d abelian duality. The abelian duality on 4d is expressed by F® v = ^t^p^F^ where 
Fpv = d^Ay — dyA^ and F® v — d^A^ 1 — d v A^, where A^,A® are the original and the dual gauge potential. 
Using dimensional reduction and splitting the duality relation into 4i and ij components, we obtain diAf = 
— \tijkFjk and F® = eijudkA^. The first relation is the well-known abelian duality relation in 3d, where we 
identify A® = a = 8 m and A4 — a D = 9 e . The monopole-instanton amplitude is naturally expressed in terms 
of 6 m — <j and e = A9. 
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In the classical background, the long distance dynamics is abelian and the sum over monopole- 
instanton amplitudes maps, via Poisson resummation, to a sum over the electric charges of 
the dyons. 

From this point of view, we find (3.21) quite intriguing. This relation makes perfect sense 
in M = 2 Seiberg-Witten theory [23]. In that context, on M 4 , the theory has a classical moduli 
space along which gauge symmetry is reduced to U{1) due to adjoint Higgsing by the vev 
($) = vT 3 . In the semi-classical domain of the N = 2 theory, the Poisson duality relation is 
reliable. Moreover, in N = 2 theory, the combination J g2 v{E + iB) ■ = Z( nmi „ c ) appears 
as the central charge in the supersymmetry algebra. 

In pure Yang-Mills theory, the combination of electric and magnetic charges J s2 (E + iB) ■ 
cE (without the vev) appears naturally. However, we also know that there is no apparent 
adjoint Higgsing on M 4 and it seems hard to make sense out of the Poisson duality relation 
in that case. One can speculate that there may be a connection between Poisson duality 
and the maximal abelian gauge proposal by 't Hooft [25], which is based on using a specific 
gauge condition applied to a composite operator to define an adjoint Higgs field, and the 
corresponding abelianized SU(N) — > U(1) N ~ 1 theory. In that case, of course, monopoles (as 
well as dyons) appear. However, the theory is strongly coupled, and the abelian dynamics is 
not well understood. This is a direction worthy of further pursuit. 

4. Comments on the literature 

In this Section, we will attempt to clarify the relation between our work and previous ana- 
lytical approaches to the deconfinement transition in the literature. 

The question whether one can study the deconfinement phase transition at weak coupling 
was considered by Aharony et al. in [26]. These authors found a realization of this idea in 
a finite spatial volume § 3 x S 1 , but rejected the possibility that an example can be found 
for an infinite spatial volume, in particular on M 3 xS 1 . Ref. [26] studied large- N (strictly, 
N = oo) pure Yang-Mills theory on small § 3 x S 1 by integrating out the perturbatively 
weakly coupled modes and thus mapping the field theory to a matrix model. This approach 
pushes the deconfinement transition to the weak coupling regime. However, since they study 
gauge theory on a small sphere § 3 , approaching the thermodynamic limit requires taking the 
infinite- N limit. In particular, the approach of [26] does not apply to gauge theories of finite 
rank. 

A way around the obstacle of analytical "intractability" of deconfinement in an infinite 
spatial volume was found in Refs. [27, 28] by compactifying the gauge theory on R 2 x x 
where is a spatial circle and §>^ is the thermal circle. The size of the spatial circle provides 
a tunable control parameter, the counter-part of S 3 in the approach of [26]. The advantage 
of this formalism is that the small volume theory is still a field theory in an infinite spatial 
volume (as opposed to a matrix model), although it is M 2 and not R 3 . The approach of [27,28] 
also works for finite rank gauge theories, mapping the deconfinement transition to a phase 
transition in two-dimensional spin systems; the relevant spin systems are the affme XY-spin 
models with symmetry breaking perturbations. 
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In the present paper, we gave a reliable semi-classical analysis of the center-symmetry 
changing phase transition on I 3 x S 1 . The use of semi-classical methods in connection with 
the center-symmetry changing transition was previously investigated in a series of papers by 
Diakonov and collaborators [20,24,29,30]: 

1) The first of these papers, a review published in 2002 [29], suggests that center- 
symmetry can be stabilized by monopole-instantons (dyons, in the language of Diakonov 
et al.). Diakonov obtains the term V ~ cosh p-(vr — AO) in the monopole-instanton induced 
potential, see (2.33) and (3.6), and observes that this contribution has a minimum at the 
center symmetric point A9 = ir. However, this conclusion is based on neglecting the the 
cos a term which arises from the coupling to the dual photon. Indeed, we found that in the 
semi-classical case m< A monopole-instantons favor center-symmetry breaking. 

2) M = 1 SYM theory on M 3 x S 1 is presented as an example for center stabilization 
by monopole-instantons in Refs. [20,24,29]. This interpretation is common in the literature 
but, as we showed above, it is physically not correct. Monopole-instantons in N = 1 SYM 
have fermion zero modes, and they do not generate a bosonic potential for the Wilson line. 
Rather, they generate a fermion bilinear which determines the superpotential. The bosonic 
potential can be found via supersymmetry using V ~ Ii5£~| 2 ) as i n Section 2.2, or using 
the BZJ-prescription, as in Section 2.3. In either case the conclusion is that the physical 
mechanism generating the potential for the Wilson line is related to topological molecules. 

3) A quantitative theory of a confining ensemble of monopole-instantons was proposed 
in [24] . The paper argues that one can analyze the vacuum of Yang-Mills theory by semi- 
classical means, and that, furthermore, the only topological configurations that contribute to 
the non-perturbative potential for the gauge holonomy satisfy: i) the self-duality condition, 
and ii) a (magnetic) charge neutrality constraint. This means that, contrary to item 1), 
instanton-monopoles [.Mi], [M2] and [Mi], [M2] are excluded. The self-dual neutral objects 
are KvBLL-calorons [8,9], 4d BPST instantons at finite temperature and non-zero holonomy, 
and multi-calorons. Calorons can be viewed as topological molecules of the form [.M1.M2] and 
[.MiA^]- We observed in (3.12) that at the classical level there is no coupling of instantons 
to the holonomy. The claim in Ref. [24] is that a potential for the Polyakov line can arise 
from the collective coordinate measure. This contribution is formally of higher order in the 
coupling constant. 

In the pure gauge theory the critical temperature for deconfinement is of order A, and it is 
clear that there cannot be a systematic semi-classical theory of the transition. The continuity 
argument outlined in Figure 1 provides a less ambitious program: We connect the strongly 
coupled center-symmetry changing phase transition in Yang-Mills theory to a semi-classical 
phase transition in mass-perturbed super- Yang-Mills theory on M 3 x S . The semi-classical 
calculation at small m provides two important lessons: a.) Both monopole-instantons (self- 
dual and magnetically charged objects 8 ) and neutral bions (non-self dual and magnetically 



It is clear that there cannot be a general argument that rules out contributions from magnetically charged 
objects. In particular, Polyakov's solution to the Yang-Mills adjoint Higgs system on R 3 [31] and the solution 
of deformed Yang-Mills theory [32] on R 3 x S 1 map the gauge theory partition function to a grand canonical 
ensemble of magnetic charges. 
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neutral objects like [.M1.M1]) contribute to the potential for the gauge holonomy and b.) 
these are the leading contributions to the potential in the controllable small-L, m regime. 
Thus, we believe that our results contradict the assumptions in [24]. 

Finally, we also note that detailed phenomenological studies of the effective lagrangian 
for the Polyakov line in pure Yang-Mills theory were carried out by Pisarski and collabora- 
tors [33-36] . Center-stabilizing double-trace operators are considered in Ref. [32] to address a 
semi-classical mechanism of confinement and large- N volume independence, and in Ref. [37] 
to study phases with partial center symmetry breaking. The present study provides a micro- 
scopic explanation of the origin of the center-stabilizing double-trace terms in these studies. 

5. Conclusions and outlook 

In this paper, we argued that the center-symmetry changing phase transition in thermal Yang 
Mills theory is continuously connected to the phase transition in softly broken N = 1 theory on 
R 3 x 8 . We showed that for small values of the adjoint fermion mass m the critical scale L c is 
analytically calculable. We also provided theoretical evidence that the same mechanism that 
drives the phase transition at small m, the competition between center-stabilizing topological 
molecules and center-destabilizing monopole-instantons, also exists in pure Yang Mills theory. 

There are a number of directions worthy of further pursuit, including both numerical and 
analytical studies: 

1. The phase diagram in Figure 1 can be explored with lattice methods that are available 
today. Furthermore, lattice simulations can be used to study the role of neutral non-self- 
dual topological defects (such as the neutral bions), for example, by using the techniques 
of Ref. [38]. 

2. One may generalize the semi-classical study of center-symmetry changing phase transi- 
tion on M 3 x S 1 to all gauge groups, exploring the role of various topological excitations 
in the symmetry realization and the nature of the phase transition. 

Let us also make a few comments on the SU(N), N > 3 case. Here, there are TV-types of 
monopole-instantons, associated with the (N — 1) simple roots and one affine root. As in 
the SU (2) theory, in the softly broken N = 1 theory, these monopole-instantons, as well 
as perturbative one-loop fluctuations, tend to break center-symmetry. There are TV-types of 
magnetic bions and iV-types of neutral bions. The magnetic bions, for all N > 3, unlike the 
N = 2 case, lead to a center-destabilizing potential. As in the SU(2) case, for all N > 3, the 
neutral bions lead to a center-stabilizing potential. For N > 3 we found a first order phase 
transition consistent with lattice gauge theory results. 

Finally, we note that the dependence of T c on the topological angle 9 can also be studied 
within our approach, resulting in a decrease of T c , in agreement with recent lattice studies 
[39] . Details of the calculations including ^-dependence and higher rank groups will be given 
elsewhere. 
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A. Supersymmetry and the non-cancelling nonzero-mode determinants 

In this work, we have derived the effective lagrangian for the Polyakov line and the dual 
photon by expanding around the supersymmetric limit of a massless Weyl fermion. In deriving 
equ. (2.14,2.33), we have used the cancellation of nonzero mode determinants in the N = 1 
supersymmetric theory. This cancellation has also been used in many other instances in the 
literature, for example in the calculation of the exact superpotential and gluino condensate 
in N = 1 supersymmetric gauge theory on M 3 x S 1 [12,40]. 

In this Appendix we discuss an important issue that has not received proper attention 
within the context of N = 1 supersymmetric theories on K 3 xS 1 — the fact that the nonzero 
mode determinants in the field of a monopole do not precisely cancel. While we will argue 
that this phenomenon has no significant effect on our main result, we include this discussion 
for completeness, as it fills a gap in the literature on JV = 1 theories on M 3 x S . 

The point is that in the supersymmetric theory, the determinants of nonzero mode fluctu- 
ations around the BPS or KK monopole-instantons on M 3 x S 1 do not precisely cancel, despite 
the fact that the solutions preserve one-half of the supersymmetry. This noncancellation oc- 
curs essentially because of the slow fall-off of the monopole-instanton background at infinity 
in R 3 [42,43]. The nonvanishing of the nonzero mode determinants around supersymmetric 
monopole-instanton backgrounds in N = 1 theories on M 3 x S 1 is the N = 1 counterpart of 
the mass and central-charge renormalization 9 of BPS monopoles in 4d iV = 2 supersymmetric 
theories. As we will see, the non-cancellation of the nonzero mode determinants is perfectly 
consistent and is, in fact, required by supersymmetry and holomorphy. 

In the context of purely 3d theories with extended supersymmetry, obtained by reducing 
N = 2 4d theories, this non-cancellation has been known since [41]. However, it was not 
addressed in the original calculation of nonperturbative effects in iV = 1 theories on M 3 x S 1 
[12,40]. Only recently, a relevant calculation on M 3 x S 1 appeared in the literature [45], 
in the context of theories with N = 2 supersymmetry. While N = 2 theories share many 
features with the N = 1 theory of interest to us, the most important difference is that the 
branch of moduli space considered in [45] is the 4d Coulomb branch, where the abelianization 
of the gauge group is due to a nonzero expectation value \a\ of the adjoint Higgs matter 
supermultiplet. In particular, |a| 3> Ajv =2 was required for consistency of the calculation 

9 See [44] for a recent review and references. 
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of [45] (while an arbitrary \a\L was allowed). In contrast, we are interested in the N = 1 
theory where the adjoint Higgs field is absent and the abelianization is due, instead, to a 
Wilson line expectation value. Thus, the result of [45] is not directly applicable to the case 
of interest to us. 

From the discussion in the previous paragraph, it is clear that a calculation relevant for 
N = 1 theories on M 3 x S 1 is not present in the literature. In what follows, we shall perform 
this calculation. The novel ingredient that we will use is the Nye-Singer index theorem on 
M 3 x S 1 [10] in backgrounds with nontrivial holonomy, in the form studied by two of us (EP 
and MU) in [11]. We will also benefit from insight gained from refs. [41,45,49]. 

The main object of interest turns out to be the "index" I(M 2 ), defined 10 as: 

M 2 M 2 
/(M 2 ) =tr A ~ tr (A.l) 

y J A_ + M 2 A + + M 2 y ' 

where A_ = D^D = -D^D* - ^a^F^ and A + = DD^ = -D^, where D is the Weyl 
operator and — the covariant derivative in the monopole-instanton background. In a self- 
dual monopole-instanton self-dual background, it is the operator A_ that has zero modes; 
see, e.g., [46] for a review of this notation. To establish a relation between I(M 2 ) and the 
nonzero mode determinants, we note the identity: 

A 2 

dM2 I(M 2 ) = trln A A + + /f - trln *~ ±£ = Indet A + + A ~ + * 



M 2 A+ + A 2 A_ + A 2 A+ + A 2 A_ + A 2 ' 

(A.2) 

where Apy is the Pauli-Villars mass and fi is an auxiliary parameter which will be eventually 
taken to zero. We can now use (A.2) to define: 



R = lim 



( a pv 2 \ 



3 

det A + det A_ + A 2 PV \ 1 



det' A_ det A + + A? 



PV 



(A.3) 



2 



Here, "det'" is the determinant with the zero modes omitted, i.e., det'A_ = lim det A a +At 

using the fact that A_ has two zero modes. The quantity R of (A.3) is, as is made clear from 
the second identity, equal to the ratio of nonzero modes' determinants around monopole- 
instanton backgrounds appearing in the N = 1 theory. Notice also that R is equal to the 
fluctuation determinant [45] for the N = 2 theory to the power of 3/2. The various con- 
tributions making up the ratio in (A.3) are as follows: the adjoint fermion contribution is 
(det' A_ det A + )*, while, in a 4d background Lorentz gauge, the gauge field determinant is 
(det' A_) _1 and the ghost determinant — (det A + )2, see [46] for details. 

10 We use quotation marks, since the value of the function 7(0) is, indeed, the index of the adjoint Dirac 
operator in a monopole-instanton background, but the quantity I(M 2 ) itself depends on M 2 [42]. 
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Our main goal here is to compute the ratio of determinants R (A. 3). We will use the 
expression for I(M 2 ) from [11]: 



I(M 2 ) = h+h{M 



2np . 2np 
2 _ *±j v . ST^ \ L ' L 




h = • (A.4) 



7T 

r2 i 



The definition of I2(M ) is evident; see also (A. 6) below. The relation between v, the ex- 
pectation value of A4, and A8, the angular distance between the eigenvalues of the Polyakov 
loop, is: 

AO = Lv G (0,2vr) . (A.5) 

We note that = 4Q, where Q is the topological charge of the solution, which equals 1/2 
at the center symmetric point v = |J. 

If L is set to zero and the p-sums restricted to p = 0, this is exactly the quantity J(M 2 ) ]R 3 , 
found in [42], and yielding the well-known result 1(0) = 2 for the adjoint representation. Note 
that the function I(M 2 ) K 3 played a role in both [41], where it was used directly, and [45], 
where the related difference between the density of states of A_ and A + , see [43], was used. 
The novelty on M 3 x S 1 is that a sum over Kaluza-Klein modes and a contribution of the 
topological charge to the index (the non integer I± term in (A.4)) appear. The first term in 
(A.4) is the bulk contribution to the index and the second — the surface term. The surface 
term is given by the Kaluza-Klein sums in (A.4), defined using zeta- function regularization; 
in fact, the second term in (A.4), 

J 2 (M 2 )= r i r , (A.6) 

as shown in detail in [11], leads to: 

vL 

-27T- 



m- t (^( 2 5? + „)-(^-„)) = - 2 ^-2 [ -g J+ 2 L ; 

p= — OO \ \ / r 



(A.7) 



where [xj is the largest integer smaller than x. Thus, the topological charge contribution to 
the index, I±, is canceled by the non integer contribution from the KK sum, -^(O), and the 
index equals 2 for < vL < 2tt. 

A.l Calculating the ratio of nonzero mode determinants 

In the previous Section, we showed that ^(0) = 2. Here, we use the expression for l2(M 2 ) 
to compute the ratio of non-zero mode determinants R. We begin with the expression for 

R given in (A. 3), multiplied by e~ So = e » 2 ( a fv) ; where So is the classical action of a BPS 
monopole: 



"So^ = e So lim 



( ,2 AP r V '2{M 2 ) ,., 2 \ 

V / 
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V(A PV )~ 



log A j 



lim 



ji A e 



-hlo gf i 2 L 2 



- I 



dM 2 



(A. 



e 9 2 (Vi) lim 



-I 1 \og^L 2 + A J V l2 ^dM 2 

/i e 



V 



e - 9 2 d/i) i? 2 



/ 



where i?2 is implicitly defined in the last line above. In other words, the UV divergent 
contribution to R serves the purpose to renormalize the coupling from h-pv to the scale 
Any other UV divergence at one loop would be in need of a counterterm and there is not 
another one at one loop. 

Next, we consider \0gR2 (omitting the explicit mention of the fi — > limit to be taken 
at the end and noticing that the upper limit of the integral can be taken to infinity): 



3-Zi 3 
log R 2 = log/i 3 - — log/iL + - 



dM 



2-rrp 



+ V 



;2-P +u )2 +M 2| 



(v 



(A.9) 



Next, noting that J 



00 dM 



/•» MV A 2 +M 2 \A\' 



arcsinh^ = In + 0(fi 2 ), we find: 



logi? 2 = log/u 3 - 3 ^log/xL+ |^ ( sign( 2 ^ + v) lo{ 



2 \ 2 -^ + v\ 



3 3/i _ 3 
log A* - — log/iL + - 



^ sign(p + 



2tt' 



(v 



log 



(0 

47T 



E. , vL.. . vL 
sign(p+ — )log|p 



2vr 1 



(v — > —v) 



(A.10) 



Before continuing, we recognize from (A. 7) that one of the KK sums appearing in (A. 10), 
'j — > —v) = 2 — Ji, giving rise, after substitution in (A.9), and remembering 



E p sign(p+ ^ 
1 



J, _ 2Lv 



log R 2 



3h 
2 



log/^ 3 - ^ log/xL + £(2 - /1) log ^- + 



4vr 3 
+ 2 



sign(p + — ) log |p + — I - (v 



2tt' 



, 4vr Lv, 3 
3 log — - 3 — log4vr+ - 

L 7T Z 



2^sign(p+ — )log|p + 



2vr 1 



(A.ll) 
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Now we can deal with the remaining KK sum, by writing it as an s-derivative of a function, 
evaluated at s = 0: 



log i?2 



, 4vr Lv, 3d 
3 log — -3 — log4vr- - — 
L 7r 2 as 



sign(p+ ^ 

p 1^ ^ 2tt I 



(u — > —v) 



Next, we define the function: 



\^ sign(p + a) 

His,a)= > — : : — 

^ \p + a\ s 



s->0 



(A.12) 



(A.13) 



and rewrite (A.12) as: 

\ogR 2 = 31og— -3— log4vr- -km — [H s, — - i? (s,- — 

L IT 2 s^O as \ Z7T 2"7T 

Now, for 1 > a > we have: 



H(s,a) = Y— i —-Y ] 1 



C(s,a) - C(s, 1 - a) 



p>0 |J 1 p>0 

where £(s, a) = Z~2 P >o \P + a l~ S ^ s ^ e incomplete zeta function, and similar for: 
H ^ "«) = E | p+ i_ a |. " E = C( S , 1 - a) - C( S , «) 



p>0 



p>0 



Thus, i?2 is given by: 



log i? 2 = 3 log — - 3— log 4vr + 3 C (0, l-—)-C(0 

L TT \ 



2vr' 



2vr 



(A.14) 



(A.15) 



(A.16) 



(A.17) 



which, upon plugging into (A. 8), and using £'(0, x) = lnT(x) — | log27r yields for e 50 times 
the ratio of determinants R in the BPS monopole background: 



~ S °R 



BPS 



4tt 

T 



4^\ 3 e -^^|y+3iogr(i-^)-3iogr(^) 



+31og47r]+3(C'(0,l-|2)-C'(0,|^)) 



For the ratio of determinants around the KK monopole, we replace 

3 



vL . i vL 

2it ~^ 1 



e~ Sa R 



KK 



-^4y + ^^wy- 31 °g r ( 1 -^)+ 31 °s r (t 



L J 



2tv , 



(A.18) 
^r, yielding: 

(A.19) 
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We can now study the expansion of the T functions from (A. 18) and (A. 19) near the 
center symmetric points: 



31og £ w 1 = ~ 3Mm " + ° (< ^ _ 1)3) 

w 5.8903 - lj , (A.20) 

where we used the value ■0(1/2) = r'(l/2)/r(l/2) ~ —1.96 of the digamma function. Thus, 
we can rewrite (A. 18) as: 



e~ s °R 



BPS 



(A.21) 



as well as a similar expression for the KK monopole (A. 19). Thus, near the center symmetric 
point, we have that there is a small shift of the scale of the coupling constant away from 
4tt/L (using —1.96 ~ — log 7.12, to a scale a few times lower than 1/R = 2ir/L) and that the 
exponential in e~ So R is a linear function of the deviation from the center symmetric vacuum, 
^ — 1, up to cubic terms. The behavior of the ratio of one- loop determinants is illustrated 
for general vL in Figure 2. 

It is of interest to also consider the neighborhood of vL = 0, in particular, if we wish to 
take the 3d limit and compare with previous studies as a check on our calculation. In the 
3d limit, we take L — > 0, but keep g\ = Ljg\ fixed; as before, v is an arbitrary position on 
the moduli space, which becomes noncompact in the 3d limit. The contribution of the KK 
monopole vanishes in the 3d limit, as its action becomes infinite. Hence, we concentrate on 
the BPS monopole amplitude. In this limit, we have that: 



3 log 



so we can write for (A. 18): 



= -31og— + 3 7S — + Y), (A.22) 

LVKO LV IT TT 



So 



R 



/47rV "57^-3 log If +3 7e ^+0((^)3) 

o 9 { L ) 



BPS 



4ttv 



-> {2v) 3 e 4 , with L -> 0, fixed g\ = — j— . (A.23) 

9i\ l / 

We note that in this limit, our calculation exactly reproduces the 3d result of [41]. To facilitate 
the comparison, we need to take into account the facts that v = my/ in our normalization 
and that our calculation is in the 4-supercharges theory, not the 8-supercharges one, which 
changes the overall power of the dimensional pre factor from 2, as in eqn. (26) from [41], to 
3 as in our (A.23). It is clear from the above that the noncancellation of the determinants is 
required in order to have a smooth 3d limit. 
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Figure 2: The behavior of the ratio of nonzero mode determinants around a BPS monopole-instanton, 
eqn. (A. 18). We define fix) = 3(logr(l — x) — logr(x)), where x — vL/2ir. The singularities of f(x) 
near the gauge-symmetry restoration points (x = 0, 1) arc logarithmic and arc, in fact, required in 
order to take the proper 3d limit, see discussion around eqn. (A. 23)). Notice that at the center- 
symmetric point the nonzero modes determinants cancel exactly, /(1/2) = 0, while the slope there is 
-3^(1/2) ~ 5.89, see eqn. (A.20). 



A. 2 Interpretation 

The main question is, what is the interpretation of the cubic and higher terms in f(x) away 
from the center symmetric point? It appears at first sight that they do not belong in the 
superpotential, as they are not holomorphic functions of vL+ia. Further, any nonlinear terms 
in this variable appearing in the exponent would be in conflict with symmetry arguments. 

We note that these questions arise already in the 3d supersymmetric theories and have 
been addressed in some cases. Notably, in the 8-supercharges 3d theory, the BPS monopole- 
instanton induced terms give rise to the four-fermi interaction [41]: 



4:7TV I - _ 

.2 ,2 ^ +WT 



Since the 8-supercharge theory admits no superpotential, this term, including the f-dependent 
prefactor should be interpreted as arising from a component expansion of a supersymmetric 
sigma model with a hyper-Kabler metric. In fact, as shown in [41], the above four-fermi term 
is fully consistent with the semiclassical expansion of the Atiyah-Hitchin metric. 

On the other hand, in the 4-supercharges 3d theory, we have, from (A. 23), after incor- 
porating bosonic and fermionic zero modes, the following form for the monopole-instanton 
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induced fermion vertex: 



(si) 3 



aV 



»:1 



-+ jo- 



lt is natural to ask how this term is absorbed in an e x superpotential, where X = — icr, 

93 

as argued long time ago [48]. The pre-exponential factor there and in subsequent work was 
assumed to be a modulus-independent constant and has not been calculated before, to the 
best of our knowledge. Thus the question regarding the incorporation of the nonholomorphic 
~ v 3 prefactor in a supersymmetric effective lagrangian applies to the 3d 4-supercharge theory 
as well. 

We now proceed to address this issue by first writing down the fermion vertices that 
accompany the BPS (A. 18) and KK (A. 19) monopole-instantons, performing the same steps 
as in [12]. We skip the details (of including the collective coordinates measure, etc., as they 
are identical to [12]) and note that the only difference is that we replace the Pauli-Villars 
scale appearing in the fermion measure there with the ratio of determinants that we calculated 
above. The fermion bilinear terms — the kinetic term and the ones that arise due to the BPS 
and KK monopole-instantons — are (all couplings below are taken at the scale -#): 



L 



ferm 



L 



Xa m d m X 



(A.24) 



+ 



2 2 4 tt 2 L /47TY 
9 2 \L) 



• j2 ' J V 27T I I g 3 2 9 2 27V 



XX 



where 



rfi 



T(x) 



(A.25) 



was already defined in the caption of Figure 2. For future interpretation, it is also useful to 
recall the original, not resummed, expression for f(x). This follows from (A. 11) and, as given 
below, is valid for < x < 1: 



f(x) = 3 In x + 3 ln(p + x) — ln(p 

p>0 



(A.26) 



Next, we rewrite (A.24), using b = p-^L and expanding around the center-symmetric value 
b = ^p- + b' (in other words, b' is the fluctuation around the center symmetric vev): 



L 



ferm 



L rfXa m d m X + Irr^wr d m b'd m b' + ... 



2 16vr 2 L 



(A.27) 



+ - 



2 10 7T 5 



e » 



XX , 



where we included the kinetic term for b' to leading order (subleading contributions to the 
b' kinetic term, to be elaborated below, are denoted by dots). The main question we want 
to address is how (A.27) can be incorporated in a supersymmetric effective lagrangian, given 
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the extra, seemingly nonholomorphic, dependence of the fermion 't Hooft vertex on b' — ia, 
arising from the noncanceling determinants around the BPS and KK monopole-instantons. 

For the reader not so interested in the details, the summary of the discussion of the 
following section is that after properly performing the supersymmetric "photon-dual photon" 
duality transformation, the dimensionless chiral superfield B dual to the abelian vector mul- 
tiplet (in 3d, equivalent to a real linear multiplet), describing and the 3d photon, is such 
that its lowest component is: 

B\ = b'-f(^ + £^b^-ia. (A.28) 

In terms of the chiral superfield (A.28), the superpotential is holomorphic and is of the well- 
known affine-Toda form, W ~ e~ B + e B . In the next Section, we discuss the details of 
the duality transformation and the interpretation of the function / in (A.28) as encoding 
perturbative corrections to the moduli space metric on M 3 x S . We recall again that / 
is nonzero only away from the center-symmetric point, thus the extremal point of W(B), 
B = 0, still corresponds to the center-symmetry preserving vacuum b' = 0, in accord with 
the conclusions from the earlier work [12,40]. 

A. 3 Linear-chiral superfield duality 

The dimensional reduction of an abelian 4d vector multiplet to 3d is described as a real linear 
multiplet, W, defined as: 11 

W=-a 36a D & D a V, (A.29) 

where V is the usual real vector superfield in 4d (dimensional reduction in (A.29) as written 
is performed along the x 3 direction and, hence, the lowest component of (A.29) is W\ = A3). 
The field W obeys the linear multiplet relations D 2 W = D 2 W = and is invariant under the 
usual (3d) supergauge transformations V — > V — i(A — At). The lowest component of W, as 
already indicated, is the real scalar in the 3d vector mulitiplet (i.e., is related by rescaling to 
our fields b or b'), while the other terms in its superspace expansion involve the fermions and 
the i7(l) gauge field strength. As defined above, the dimension of W is unity. The minimal 
kinetic term for the vector multiplet is given by the first term in the D-term action given 
below: 

J d 4 0d 3 x (~^W 2 - F(W)j > ( A - 30 ) 

while the function F(W) incorporates possible nonlinear corrections to the kinetic terms. 
The action (A. 30) is the most general two-derivative one involving the linear multiplet (note 
that W can not appear in integrals over half-superspace). The coupling e 2 in (A. 30) has unit 
mass dimension and denotes the 3d gauge coupling; clearly, e 2 or another mass scale must be 
present in F(W). The kinetic term of the lowest component of W, which we denote by w in 
this Section (W\ = w), hoping not to cause confusion, reads: 

/ d3x F "^ d ^ dmyj ■ ( A - 31 ) 

n Our notation, including the supercovariant derivatives and V, is that of Wess and Bagger [47]. 
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For the theory to make sense, the function \ +F"(w) (the "tau-parameter," or moduli space 
metric), which also determines the effective U(l) coupling along the Coulomb branch, has to 
be positive. 

Linear-chiral duality is performed via a Legendre transformation as follows: 

J d A 9d 3 x |-^ 2 - F(W) + (fit + B)w) . (A.32) 

The dimensionless chiral superfield B is introduced as a Lagrange multiplier and W is now 
regarded as an unconstrained real superfield. Variation with respect to B and B\ taking into 
account the chirality constraint D a B^ = D&B = 0, enforces the linear multiplet conditions 
D 2 W = D 2 W = on the real superfield W. 

On the other hand, a variation of (A.32) with respect to W gives: 12 

B + B ] = \w + F'(W) . (A.33) 
e z 

The equation of motion (A.33) can be inverted, e.g. perturbatively, to give W = e 2 (B + 
B*) + . . ., and upon plugging the solution of (A.33) in (A.32), one is left with a dual theory 
of chiral superfields given by: 

J d A 6d 3 x K(B + £t) = J d 4 9d 3 x (j(B + B^f + = J d A 9d 3 x [e 2 B^B + . 

(A.34) 

where the dots denote higher order terms determined by the form of F(W). For a recent 
reference, see, for example, [49]. The point of the discussion above is that the relation 
between the real part of the lowest component of the dual chiral superfield B and the lowest 
component of the linear superfield W is simple — i.e., linear — only if the function F(W) is 
quadratic, i.e., if there are no nontrivial corrections to the kinetic terms of the b' ~ w field. 
In general, this relation can be complicated. 

In particular, in our case of interest, the duality transformation on M 3 x S 1 is performed 
after Kaluza-Klein modes are integrated out. Integrating out the KK modes generates non- 
trivial corrections to the kinetic terms of A4 and 3d photon. As we will now see, these 
perturbative corrections are precisely related to the function / (A. 26) representing the ratio 
of nonzero modes' fermionic and bosonic determinants around the BPS and KK monopole 
instantons. To see this, let us now apply the duality transformation to our theory. As already 
noted, eqn. (A.33) determines the lowest component of X in terms of w(x) and the dual 
photon a(x): 

B\ = ^w + ^F'(w)-ia , (A.35) 
which, upon comparison to (A. 28) immediately implies that: 

b'(x) = ^w(x) and F>(w) = -2/ Q + j^w) ■ (A.36) 

12 Notice that this relation also implies that the fermion components of B are slightly different from the 
fermionic components of W; we will ignore this in what follows. 
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Then, we can use (A. 36) to compare the leading order kinetic term of w, eqn. (A. 31), with 
that of 6', eqn. (A. 27), and find: 

.2 9 2 



- - UM ' < A ' 37 > 
We can now use the explicit form (A. 26) of / and (A. 37) to find that the corrections to 
the kinetic term for w (A. 31) are determined by the effective "tau parameter": 



e z 
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x p — x p + X 



(A.38) 
(A.39) 



In the last line, we used eqn. (A. 26) for f(x) and denoted f'(x) = df(x)/dx. 
In order to interpret the corrections, we rewrite (A.38) as: 
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(A.40) 



It is most straightforward to first consider the purely-3d limit L — > 0, Jj-fixed. In this 
case, the sum over p > in (A.40) drops out and only the first term in the square brackets 
survives. Recall that g 2 is taken at the scale ^7, which is now the 3d UV cutoff scale (the 



bare 3d coupling is p 



I. The coupling (A.40) can be rewritten, using x 
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(A.41) 



The shift of i n (A.41) represents exactly the one-loop shift of the bare 3d gauge 
coupling due to integrating out the 3d heavy vector multiplet — the heavy VF-boson and 
superpartner, of mass ~ v, along the Coulomb branch parameterized by v in an SU(2) theory. 
This one- loop shift was calculated in [50], see eqn. (83) there. 13 Note that we obtained the 
same result as [50] by a rather roundabout way — by calculating the ratio of determinants of 
nonzero modes in the monopole-instanton background and then demanding consistency of 
the result with holomorphy of the superpotential. While this may appear miraculous, the 
matching of the results had to be true, by the power of supersymmetry. 

Going back to finite-L, i.e. to M 3 x S 1 , the interpretation of the other terms appearing 
in (A.40) is now clear. The sum over p > represents the effect of the Kaluza-Klein partners 
of the heavy vector multiplets on the unbroken- U(l) gauge coupling; recall that at every KK 
level with p > 0, there are two heavy vector multiplets of mass 2 j^±v. We are not aware of an 
explicit perturbative calculation of this in a compactified TV = 1 theory however, similar to 

13 The factor of 64-7T 2 represents an overall normalization factor; note that ref. [50] calculated precisely the 
"running" of the 3d coupling from the bare value 4r to the value at scales below v, with a result exactly as 
given in the brackets in (A.41). 
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the purely-3d four supercharge case and to the N = 2 4d theory, and given the suggestive 
form of the KK sum in (A. 40), it is natural to conjecture that, by supersymmetry, the 
nonvanishing ratio of nonzero modes determinants around the monopole-instanton solutions 
precisely encodes the one- loop corrections to the U(l) gauge coupling. 

In our discussion in the main text, we have ignored the non-canceling nonzero mode 
determinants. This is because their net effect is to modify the Kahler potential for the 
chiral superfield B, see eqn. (A. 35), by small loop-suppressed terms. Our main concern is 
the competition between nonperturbative effects at weak coupling, while the perturbative 
deformation of the Kahler potential (2.10) provides only small corrections to the leading 
order semiclassical results. 



B. The monopole- induced potential in softly-broken SYM 

In this appendix, we determine the prefactor of the monopole instanton amplitude in softly 
broken J\f = 1 theory on M 3 x S . From [12], the collective coordinate measure for Mi (BPS) 
monopole-instantons in SU{2) gauge theory, is: 

d^BPS = J d\ d^ d\ e- b+ - , (B.l) 

where b = (b) + b(a), a = (a) + a (a). The three-dimensional vector a denotes the collective 
coordinate representing the center of the monopole, <fi is the angular U(l) collective coordinate, 
£ a are the two Grassmann fermion zero mode coordinates, and the prefactor represents the 
product of all collective coordinate Jacobians. 

Now, we add the fermion mass term AL m = ptr[AA], see (2.32), to the Lagrangian 
of SYM. In order to obtain the contribution of, for example, the Mi (or BPS) monopole- 
instanton to the potential for b and a, we saturate the integral over the fermion zero modes 
by a single mass-term insertion: 

J dfi B pse- ASm « -Mh 1 J ffae-b+ia mL J ^ ^ _ _ ^ (R2) 

where AS m = f d 4 x AL m and X°(x — a) is the fermion zero mode of the monopole located at a. 
Next we note that J d 2 ^ d?x tr[A°A°] is exactly the fermion zero-mode Jacobian calculated 

in [41], which is equal to 25^ / ' 1 = 2b. Collecting all factors, we conclude that the contribution 
of a Mi monopole-instanton to the potential is: 

^ V M, = 2mM PV L })e -b + ia (R3) 

Proceeding similarly, we find that the Mi (or KK) monopole-instanton contributes, using its 
action S^ 2 = ^ — b instead: 

AV M 2 = *nM*pyL W _ A (R4) 

g 2 V 9 J 



14 See Section 4.1 in [51] for an expression, similar to (A. 40), for the one-loop correction to the moduli space 
metric of the N = 2 theory on R 3 x S 1 , valid also on the Coulomb branch with only a Wilson line turned on. 



To obtain the total 0(m) contribution to the scalar potential, we now sum over the contri- 
butions of the two monopole-instantons and their complex conjugates: 



Finally, we use the relations (2.16), to cast V mon _ into the form given in (2.33). 
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